Coherent radiation from neutral molecules moving above a grating 
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We predict and study the quantum-electrodynamical effect of parametric self-induced excitation of 
a molecule moving above the dielectric or conducting medium with periodic grating. In this case the 
radiation reaction force modulates the molecular transition frequency which results in a parametric 
instability of dipole oscillations even from the level of quantum or thermal fluctuations. The present 
mechanism of instability of electrically neutral molecules is different from that of the well-known 
Smith-Purcell and transition radiation in which a moving charge and its oscillating image create an 
oscillating dipole. We show that parametrically excited molecular bunches can produce an easily 
detectable coherent radiation flux of up to a microwatt. 

PACS numbers: 42.50.Gy, 42.50.Md, 33.50.-j, 32.80.Qk 
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Introduction. The presence of conducting or dielectric 
surfaces near an atom or a molecule modifies its prop- 
erties in a fundamental way, changing its radiation and 
back reaction as well as the electromagnetic vacuum fluc- 
tuations, and giving rise to a number of interesting ef- 
fects; see, e. g, § g g § | § |, 1 1, @. In particular, 
the atomic energy levels and the radiative decay rate are 
changed. A boundary can alter also the dynamics of 
atomic or molecular dipole oscillations. Here we describe 
a new effect of the latter type. 

The calculations of the level shifts for a molecule 
near a perfectly conducting wall can be traced back to 
the papers Ji , [Tl| , fl2) . According to subsequent works 
p3| , fl4| , |l5| , jl| (and in accordance with fluctuation- 
dissipation theorem), this effect is due to the modifi- 
cation of both electromagnetic vacuum fluctuations and 
radiation reaction. Note, however, that in the near zone 
R <C Ao/^/eT the influence of a boundary between two 
media with dielectric constants e± and £2 on a radiat- 
ing dipole is dominated by the radiation reaction force, 
while the contribution from vacuum fluctuations of the 
electromagnetic field is negligible. This result has been 
obtained in jl5| using the ideas of separating the self- 
action and vacuum-fluctuation terms developed in [jidj; 
see || for the review. The leading term in the self-action 
force is proportional to a large factor (X /R) 3 and can 
be interpreted as a non-retarded London-van-der-Waals 
interaction with an instantaneous image dipole. Here 
Ao = 2irc/u>o is the vacuum wavelength of a given dipole 
transition with frequency ljq, e\ is the dielectric constant 
of the medium in which an atom is located. This large 
factor originates from the near field 
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created by a high-frequency image dipole p' at the posi- 
tion of a real dipole p. Here and below upper and lower 
signs correspond to the dipoles oriented perpendicular 
(_L) and parallel (|j) to the boundary. 

There is a number of reasons why the dynamics of 
molecular radiative transitions is modified in the pres- 
ence of boundaries. The effects that were identified 
and observed include (i) changes in the spectral den- 
sity of radiated modes ||, 0, ||, [| (ii) location of 
an atom in the nodes or maxima of resonant modes or in 
a non-transparent medium, e. g., in a medium with neg- 
ative dielectric constant or in a photonic band structure, 
when the transition frequency is inside the Bragg gap 
]l9| , |2(i| , pl| , (iii) phase shift of the near field of a dipole 
due to dissipation in the neighboring medium p2| , p3fl . 

Additional possibilities to affect the molecular dynam- 
ics arise when the dipoles are moving near the interface 
of two media or inside a medium. A well-known example 
is Cherenkov radiation of an oscillating dipole moving in 
an anisotropic medium under the conditions of anoma- 
lous Doppler effect which can lead to instability of the 
dipole oscillations [Q. This instability is, however, very 
difficult to realize since radiation in the directions corre- 
sponding to the normal Doppler effect usually dominates. 

In this paper we investigate a new mechanism of 
boundary-induced excitation of molecular transitions 
which is realized when a molecule is moving close to the 
dielectric or conducting medium with periodic grating. 
In this case the radiation reaction force acting on an os- 
cillating dipole moment of a given transition is mainly 
due to the time-dependent London-van-der-Waals inter- 
action experienced by the molecule moving above a grat- 
ing. It is proportional to the instantaneous value of the 
dipole moment with a proportionality coefficient being 
a periodic function of time. Its modulation frequency 
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v is equal to the velocity of a molecule divided by the 
spatial period of a grating. Such a modulated radiation 
reaction force represents a periodic perturbation of the 
transition frequency uj q (eigcnfrcqucncy of dipole oscilla- 
tions) which can drive a parametric instability under the 
conditions of Mathieu's resonance: v = 2ujq/N, where N 
is 1, 2, 3 ... In other words, the instability is due to the pe- 
riodically modulated radiation reaction force which gives 
rise to excitation of dipole oscillations even from the level 
of quantum or thermal fluctuations. 

If we have a beam of molecular dipoles, one can expect 
generation of measurable coherent radiation flux. 

We emphasize that this situation is different from the 
well-known Smith-Purcell effect |25| in which a charge 
moving above a periodic grating creates an oscillating 
dipole due to its oscillating image. It is also different 
from the transition radiation of a charge traversing a pe- 
riodic dielectric stack ptifl . In the present case, it is the 
electrically neutral atom or molecule that gets paramet- 
rically excited. 

To find the conditions for parametric instability, we 
calculate the back reaction force acting on a dipole oscil- 
lator with a high-frequency dipole moment p(t) which is 
moving with a constant velocity v in a medium with real 
dielectric constant e\ near the boundary of a medium 
with complex dielectric constant £2. Precisely, we cal- 
culate the electric field E(i) created by a high-frequency 
dipole at its position z = R, x — vt, y = 0. We assume 
that the medium 2 has a periodic grating with period 
L in the direction x along the dipole velocity, so that 
the distance R between the dipole and the boundary is 
actually a periodic function of time with period T = L/v. 

We assume that the distance R to the boundary is 
much less than the grating period, so that the radiation 
reaction field at any moment of time t is approximately 
equal to the one for the dipole at a distance z(t) above 
an infinite plane boundary. This is of course the most 
interesting limit. In the opposite case, when R^> L, the 
dipole will "feel" simultaneously many grating periods 
and the effect of grating will be averaged almost to zero. 

The problem of the radiation reaction field acting on 
a dipole above a plane boundary was considered many 
times. We will use here the results of our papers [^2|, p3| . 
According to these works, the radiation reaction field can 
be represented as a sum of a free-space term Ef ree (i) = 

2p(i) v /eT/3c 3 and a boundary-induced term E;, that can 
be very complicated. 

If the medium 2 is a perfect conductor, 62 — > ioo, then 
the contribution of a boundary is equivalent to the field 
of an image dipole p|| = — p|| or p' ± = p±, retarded by 
the time t\ — 2R^feljc: 
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This result was obtained in many works, including the 
full quantum-electrodynamical treatment of the problem; 
see Refs. [1-5,13,15,17,19,22,23]. 

To illustrate the physical mechanism of an instability 
we will assume in this paper that the medium 2 is a per- 
fect conductor and use the above expression (||) for the 
field. The case of an arbitrary medium is qualitatively 
similar and can be analyzed by using general expressions 
for the field obtained in |2^, |23| . 

Parametric instability of a harmonic dipole oscil- 
lator. This case is relevant for dipole transitions in a 
system with many quasiequidistant levels, for example, 
vibrational transitions in a molecule or Rydberg transi- 
tions in an excited atom. Such a system has essentially 
classical dynamics. The equation for free dipole oscilla- 
tions of a classical harmonic oscillator with a charge e, 
mass m, and frequency loq takes the form 



d 2 p/dt 2 + lu 2 P = (e 2 /™)(E fr cc + E b 
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In the free space, this equation describes an oscillator 
decaying with radiative rate 7 = 2e 2 cu 2 v / £T/(3toc 3 ). 

In the presence of a boundary, let us consider for defi- 
nitencss the dipole oriented perpendicular to the interface 
between two media. The result for a parallel dipole is the 
same apart from a different numerical coefficient. Then, 
expanding the electric field Ef,j_ in Eq. (||) in powers of a 
small parameter R/Xq, we obtain the following equation: 
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The factor 27 describes the well-known result that a 
perpendicular dipole above the metal surface decays with 
a rate twice that in a free space. Our main interest here 
is in the frequency shift factor which is proportional to 
1 / R? and can have much larger value than radiative de- 
cay proportional to 1/A 3 . Note that for an arbitrary 
medium 2 with complex dielectric constant £2 this fre- 
quency shift has also imaginary part cx e 2 '/ R 3 that can 
be larger than 7. This can lead to strong modification of 
the spontaneous emission rate and even reversal of the 
direction of radiative transitions [^2[ p3| . 

For an ideal conductor the frequency shift is real, and 
the only possible mechanism of instability is parametric 
excitation of dipole oscillation due to periodic resonance 
modulation of this term. This can be achieved by, e.g., 
modulating the distance R to the boundary. Such a mod- 
ulation occurs when, for example, the molecule is moving 
above a periodic grating of period L. Note that we con- 
sider here free dipole oscillations, without any external 
electromagnetic field force. 

Suppose for simplicity that the grating is sinusoidal. 
Then the distance R in the expression for the radiation 
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reaction field will be periodically modulated as 

R = i?o(l + a cos ft), 

where v = 2nv/L. The parametric resonance occurs 
when f — 2ujo/N, for an integer N = 1,2, ... . 

If the relative grating amplitude is small, a <C 1, we 
obtain from Eq. (||) the Mathieu's equation: 



p + 2jp + coq (1 + A cos ft)p = 0, 



where 
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The condition for instability on the main resonance N — 
1 reads A > Aj/luq, which can be rewritten as 
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Assuming a — 0.1, the numerical factor on the left-hand 
side of (0) is of order 10 -4 , which requires Xq/Rq > 20. 
Since i?o cannot be reasonably made much smaller than 
0.1 /urn, we need to use long- wavelength transitions, e.g., 
vibrational-rotational transitions in molecules. 

In addition, the resonance condition implies Ao = 
2cL/v. The value of L is bounded from below: L > Ro > 
0.1 //m. The velocities of neutral molecules achieved in 
supersonic jets or Laval nozzles are of the order of 1 
km/s. In this case A falls into the cm range correspond- 
ing to rotational transitions that are better described by 
a two-level model (see below). However, by accelerat- 
ing charged molecules and then neutralizing them much 
greater velocities can be achieved that allow one to em- 
ploy IR and sub-mm vibrational transitions with quasi- 
equidistant energy levels. 

Not too close to the threshold, the growth rate of dipole 
oscillations is 



uj" ~ uj q A/4 ~ 3e 2 a/(16mw i?o£i)- 
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Using the correspondence principle we can relate the clas- 
sical dipole amplitude with the matrix element d of a 
given dipole transition: e 2 /m = 2u>ocP/Tl. Then Eq. (||) 
yields the growth rate 



uj" = 3d 2 a/(8hRle 1 ). 
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For a typical dipole moment d ~ 1 Debye and Rq ~ 0.1 
/im, we get uj" ~ 5 x 10 4 s , which means that an 
excitation length v/ui" is 2 cm. Therefore, the excita- 
tion length is short enough to make the effect readily 
observable in experiment. 

Instability of a two-level system. An opposite lim- 
iting case of a molecular transition is a two-level system. 



A good example in the low frequency range is rotational 
transitions in molecules. One can obtain the modified 
Bloch equations for a two-level molecule in the vicinity 
of an interface following the usual derivation of master 
equations for a two-level system in a free space. The re- 
sult is presented below, again for a perpendicular dipole: 



p + 2 7 p + (uj 2 



d{An) jdt + 7(An - An p ) = 



d 2 uj An/(2hR 3 ))p 
2And 2 Lu h" 1 E, cxt , (10) 
-(2/nw )E cxt p. (11) 



Here An = m — ri2 is volume density of the population 
difference between two states, An p is the population dif- 
ference supported by an incoherent pumping, E cx t is an 
external (classical) field which does not contain any back- 
reaction field. 

In the absence of an external field we obtain the growth 
rate of the parametric instability, 
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which is different from Eq. ^ only by a factor An. This 
factor can be rather small for the rotational transitions at 
room temperature. This may require low temperatures 
for more efficient excitation. To reach An ~ 1/2 for 
ui{) = 10 11 s , liquid helium temperatures are required. 

The value of the growth rate can be enhanced even 
further by launching molecules closer to the surface or 
decreasing the dielectric constant e\. It is difficult to 
send molecules closer than 0.1 /im to the surface. At the 
same time, the dielectric constant can be drawn to very 
small values by creating a rarefied background plasma 
with plasma frequency close to the frequency ujq of dipole 
oscillations. For ojq — 10 11 s^ 1 the density of ionized 
gas should be around 10 12 cm~ 3 . The plasma of such 
density can be very easily produced in a gas discharge. 
Collisions in such a rarefied plasma are unimportant, and 
the residual value of ei will be defined by density close 
to the grating surface, where the double electric layer is 
formed. Experiments indicate the possibility to decrease 
El down to 10~ 4 . 

Such an enhancement of the parametric instability in a 
rarefied plasma is interesting by itself. The experimental 
realization of this effect is facilitated by the fact that 
fluctuations in E\ are not important since the quantity E\ 
in all expressions for the growth rate is actually the value 
averaged over a very large spatial scale Ai = Ao/ ' \fel- 

Possible experiments. There is a number of ways how 
this unusual dynamics of a molecule can be observed. 
One evident possibility is to detect radiation due to ex- 
cited dipole oscillations from a beam of molecules fly- 
ing above the grating. The density n in a molecular 
beam moving with supersonic velocity in a supersonic 
jet or Laval nozzle can be as high as 10 17 cm~ 3 . How- 
ever, since only the molecules moving very close to the 
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surface [L < 0.1 (ira in our numerical example) are ef- 
ficiently excited, the total amount of excited molecules 
above the grating plate of 1 cm xlO cm size will be of 
order TV ~ 10 13 . Each molecule radiates with a power 
of order W\ ~ -fhojo- At a frequency luq = 10 11 s -1 
their total spontaneous emission power will be very small: 
NW\ <~ 10~ 16 W. To increase the radiated power, one 
should pre-phase the dipole oscillations of molecules at 
the entrance to the amplification region by an exter- 
nal microwave field of intensity higher than the thermal 
noise. Then the molecules will radiate coherently in the 
bunches of length of the order of wawelength A, and the 
power radiated by one bunch will be TV 2 times the power 
radiated by a single molecule, where N <~ n\ 2 R . For 
the above values of n and A <~ 1 cm we obtain the total 
power in the microwatt range which is easily detectable. 

Another way to study the dynamics of molecules is 
to observe the change in absorption of a microwave ra- 
diation propagating through the layer of parametrically 
excited molecules. This scheme is also easily realizable. 

To increase further the radiated power and to facili- 
tate the requirements for experiments it is desirable to 
launch molecules with higher speed. This increases the 
frequency of parametric modulation and leads to greater 
power of spontaneous emission (which goes as w 4 ). Then 
one can excite vibrational transitions in molecules and 
obtain an emitter in the far-IR range where very few 
sources exist. Alternatively, for a given resonant fre- 
quency we can increase the spatial period L of the grating 
and the distance R from the interface. However, it is dif- 
ficult to accelerate neutral molecules to velocities higher 
than 1 km/s. A way to overcome this difficulty is to accel- 
erate molecular ions and then neutralize them by charge 
exchange. This can be achieved by putting a buffer gas 
cell or a thin foil before the entrance to the excitation 
region with grating interface. 

We conclude that there is one more fundamental ef- 
fect in the realm of QED, namely, the parametric self- 
excitation of dipole oscillations of a molecule moving 
close to a periodic grating. This effect can be experi- 
mentaly observed in different ways and can be even em- 
ployed for generation of coherent radiation in the far-IR 
to microwave ranges. 

Acknowledgments. We thank Marian Scully for valu- 
able discussions and appreciate the support from the 
Texas Advanced Technology Program, the Office of Naval 
research, and High Energy Lasers Program of Depart- 
ment of Defense. 



[2] R.R. Chance, A. Prock, and R. Silbey, Adv. Chem. Phys. 
37, 1 (1978); 

J.M. Wylie, J.E. Sipe, Phys. Rev. A 30, 1185 (1984); 

Ibid. A 32, 2030 (1985). 
[3] E.A. Hinds, Adv. Atom. Mol. Opt. Phys. 28, 237 (1991); 

Adv. Atom. Mol. Opt. Phys. Suppl. 2, 1 (1994). 
[4] H. Morawitz, Phys. Rev. 187, 1792 (1969); 

P.W. Milonni, P.L. Knight, Opt. Comm. 9, 119 (1973). 
[5] H.F. Arnoldus, T.F. George, Phys. Rev. A 37, 761 

(1988). 

[6] H. Khosravi, R. Loudon, Proc. Roy. Soc. Lond. A 433, 
337 (1991); Ibid. A 436, 373 (1992); 
S.M. Barnett, B. Huttner, and R. Loudon, Phys. Rev. 
Lett. 68, 3698 (1992). 
[7] A.N. Oraevsky, Usp. Fiz. Nauk 164, 415 (1994). 
[8] V.F. Cheltsov, Nuovo Cim. D 17, 17 (1995). 
[9] D. Meschede, Phys. Rep. 211, 201 (1992). 
[10] Y. Yamamoto, S. Machida, and G. Bjork, Opt. Quant. 

Electr. 24, 215 (1992). 
[11] H. Casimir, D. Polder, Phys. Rev. 73, 360 (1948). 
[12] G. Barton, Proc. Roy. Soc. Lond. 320, 251 (1970); Ibid. 

A 410, 141, 175 (1987); J. Phys. B 7, 2134 (1974). 
[13] P.W. Milonni, Phys. Rev. A 25, 1315 (1982); Phys. 

Scripta T 21, 102 (1988). 
[14] J. Dalibard, J. Dupont-Roc, and C. Cohen- Tannoudji, J. 

de Physique 43, 1617 (1982); Ibid. 45, 637 (1984). 
[15] D. Meschede, W. Jhe, and E.A. Hinds, Phys. Rev. A 41, 
1587 (1990). 

[16] E.A. Hinds, V. Sandoghdar, Phys. Rev. A 43, 398 
(1991). 

[17] A.O Barut, J.P. Dowling, Phys. Rev. A 36, 649 (1987). 
[18] G.S. Agarwal, Phys. Rev. A 12, 1475 (1975); 

G. S. Agarwal, S.D. Gupta, R.R. Puri, Fundamentals of 
cavity quantum electrodynamics, World Scientific, Singa- 
pore (1995); 

H. Morawitz, MR. Philpott, Phys. Rev. B 10, 4863 
(1974). 

[19] V.P. Bykov, Zh. Eksp. Teor. Fiz. 62, 505 (1972), Ibid. 63, 
1226 (1972) (Sov. Phys. - JETP 36, 646 (1973)); Kvant. 
Elektron. (Moscow) 1, 1557 (1974) (Sov. J. Quant. Elec- 
tron. 4, 861 (1975)). 

[20] S. John, T. Quang, Phys. Rev. A 50, 1764 (1994). 

[21] T. Suzuki, P.K.L. Yu, J. Opt. Soc. Am. B 12, 570 
(1995); M. Scalora, J.P. Dowling, M. Tossi, M.J. Bloe- 
mer, CM. Bowden, J.W. Haus, Appl. Phys. B, 60, S57 
(1995). . 

[22] A. A. Belyanin, V. V. Kocharovsky, and Vl.V. 

Kocharovsky, Laser Physics 5, 1164 (1995). 
[23] V.V. Kocharovsky, Vl.V. Kocharovsky, and A. A. 

Belyanin, Phys. Rev. Lett. 76, 3285-3288 (1996). 
[24] V.L. Ginzburg and V.M. Fain, JETP 35, 817 (1958). 
[25] S.J. Smith and E.M. Purcell, Phys. Rev. 92, 1069 (1953). 
[26] V.L. Ginzburg and I.M. Frank, JETP 16, 15 (1946). 



[1] K.H. Drexhage, in Progress in Optics, edited by E. Wolf, 
12, 163 (1974). 



